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Abstract. The relativistic quark model is presented. The quark-antiquark potential for the 
Schrodinger-like equation is constructed with the account of retardation effects and one-loop ra- 
diative corrections. It consists of the one-gluon exchange part and the confining part which is the 
mixture of the Lorentz scalar and Lorentz vector contributions. The latter contains both the Dirac 
and Pauli terms. In the v^/c^ approximation the mass spectra of heavy quarkonia (charmonium and 
bottomonium) are calculated in good agreement with experiment. In the case of heavy-light mesons 
{B and D) the light quark is treated completely relativistically and only the expansion in the in- 
verse heavy quark mass is used. The mass spectra of the ground and excited states of D, D,, B, Bs 
mesons are calculated. They exhibit some features of the so-called "level inversion". The obtained 
results are generally in accord with experimental data. Still there exist some discrepancies between 
measurements of different collaborations 



INTRODUCTION 

The heavy flavour studies lie on the frontiers of elementary particle physics. The inves- 
tigation of heavy meson mass spectra provides substantial information about the non- 
perturbative content of quantum chromodynamics (QCD) and fundamental parameters 
of the relativistic quark model. Since it is impossible to give here any kind of compre- 
hensive review of the subject we present instead the consideration based mainly on two 
papers [1,2]. They are extended to include recent experimental data and should be taken 
as an illustration of theoretical approaches. 

RELATIVISTIC QUARK MODEL 

In the quasipotential approach the meson is described by the wave function of the bound 
quark-antiquark state, which satisfies the quasipotential equation of the Schrodinger type 
in the centre-of-mass frame: 

where /jr is the relativistic reduced mass 



and b^{M) denotes the on-mass-shell relative momentum squared 



The kernel V(p,q;M) is the quasipotential operator of the quark-antiquark interaction. 
It is constructed with the help of the off-mass-shell scattering amplitude, projected onto 

the positive energy states. 

We have assumed that the effective interaction is the sum of the usual one-gluon 
exchange term and the mixture of vector and scalar linear confining potentials. The 
quasipotential is then defined by 

y(p,q;M) = M,(p)%(-p)|la,D^v(k)^Yj + Ok)r^re;^ 

+ ViM{p,q;M)^Ug{q)uQ{-q), (2) 

where is the QCD coupling constant, D^y is the gluon propagator in the Coulomb 
gauge and k = p — q; the Dirac spinor is given by 



1 



2E{p) 



with e{p) = a/p^ + wi^. The effective long-range vector vertex has the form 

where K is the nonperturbative anomalous chromomagnetic moment of quarks. Vector 
and scalar confining potentials in the nonrelativistic limit reduce to 

VcoM = (l-£)(Ar + 5), Vi,,{r)^e{Ar + B), 

reproducing V'conf('") = Konf ('") + K;onf(^) =^f + B, where 8 is the mixing coefficient. 
The potential parameters are as follows: 

A = 0.18GeV2, fi = -0.30GeV, a^(mc)=0.32, a^(m;,) = 0.22. 
The constituent quark masses are 

lUh = 4.88 GeV, m, = 0.50 GeV, - 1.55 GeV, m,,j = 0.33 GeV. 
The value of the mixing parameter £ = — 1 is fixed by matching the heavy quark 
expansion and from the consideration of charmonium radiative decays (J/'\\f ^ ^cj)- 
The anomalous chromomagnetic quark moment K = — 1 is determined from the heavy 
quark expansion and from fitting the fine splitting of heavy quarkonia ^Pj states. The 
long range chromomagnetic contribution to the potential is proportional to ( 1 + k) and 
vanishes for K = — 1 in accord with the flux tube model. 



HEAVY QUARK- ANTIQUARK POTENTIAL 

With the account of retardation effects and one loop radiative corrections the spin- 
independent potential reads as [1] 
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Here is a number of flavours, // is a renormalization scale and the subscript W means 
the Weyl ordering of operators. 

The spin-dependent part of the quark-antiquark potential for equal quark masses 
(tria = mi, = m) is given by [1] 



VsD(r) = ah-S + b 
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where L is the orbital momentum and S^^/,, S = + are the spin momenta. The total 
angular momentum is J = L + S. 



+ (l+K)^(l-8)-L 



HEAVY QUARKONIUM MASS SPECTRA 

The heavy quarkonium is similar to the positronium atom and its levels are usually 
specified by the notation n^^'^+^^Ly, where n is the radial quantum number. The results 
of calculations of heavy quarkonium mass spectra on the basis of Eqs. (1), (3), (4) are 
presented in Tables 1 and 2 [1]. 

Recently the contribution of the finite c-quark mass to the bottomonium mass spec- 
trum in one loop was considered in our model [4]. The correction to the bb static poten- 
TABLE 1. Charmonium Mass Spectrum (GeV). 



State Particle Theory PDG(2000) BES CLEO E835 

(«(25+l)Ly) II] [3] 



l'5o 


Tic 


2.979 


2.9798(18) 2.9763 2.9804 




J/^ 


3.096 


3.09687(4) 




IcO 


3.424 


3.4150(8) 3.4141 3.4154 


iM 


Xcl 


3.510 


3.51051(12) 


l^Pi 


Xc2 


3.556 


3.55618(13) 


2^So 




3.583 


3.594(5) 


2^Si 




3.686 


3.68596(9) 


I'D, 




3.798 


3.7699(25)* 






3.813 








3.815 




2'Po 


XeO 


3.854 




2^Pi 


Xcl 


3.929 




2^P2 


Zc2 


3.972 




3^ So 




3.991 








4.088 


4.040(10)* 


2^Di 




4.194 


4.159(20)* 


2^D2 




4.215 




2^D3 




4.223 





* Mixture of 5 and D states 



TABLE 2. Bottomonium Mass Spectrum (GeV). 
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State 

!2,V-I)r , 

Lj) 


JT til LiV_^iC 


Theory 

in 


|3| 




' It) 


9.400 






Y 


9.460 


9 46030^26^ 




Xbo 


9.864 


9.8599(10) 




Xbi 


9.892 


9.8927(6) 






9.912 


9.9126(5) 


2^50 




9.990 




2^51 


Y' 


10.020 


10.02326(31) 






10.151 








10.157 








10.160 






y',„ 


10.232 




23/'i 




10.253 


10.2552(5) 


23P2 




10.267 


10.2685(4) 




v!L 


10.328 




3^51 




10.355 


10.3552(5) 


2^Di 




10.441 








10.446 




22£»3 




10.450 






X/'o 


10.498 




33/^1 




10.516 




3^P2 




10.529 








10.578 




4^Si 




10.604 


10.5800(35) 



CLEO MDl 



9.46051 



9.8600 
9.8937 
9.9119 



10.0235 



tial due to mc 7^ is approximately given by [5] (flo = 5.2, = 0.5772 .. .) 

4 cc^ /"°° 
SVW = -H— NV«0"^c^)+Y£+Ei(V«omcr)], Bi{x) = / -e'' . 

" Ttf Jx t 

Averaging over solutions of Eq. (1) with the Cornell potential yields the following 

State 15 IP 25 \D IP 35 

bottomonium mass shifts: , . 

(5V),MeV -12 -9.3 -8.7 -7.6 -7.5 -7.2 

These shifts are within the estimates of theoretical uncertainties (~ 10 MeV) and par- 
tially could be adsorbed either in the value of the constituent b quark mass or in the 
constant term B of the confining potential. 



MASS SPECTRA OF B AND D MESONS 



The B and D mesons are alike the hydrogen atom with the heavy quark Q = b,c near the 
centre-of-mass and the light quark q — u,d,s orbiting around it. Thus it is appropriate 
to make the expansion in the inverse heavy quark mass. In the limit of infinitely heavy 
quark (niQ — > oo) its mass and spin decouple and as a result heavy quark symmetry arises. 
In this limit the meson angular momentum is a sum of the orbital momentum L and the 
light quark spin i.e. j = L + S^. The 1/mg correction to the Qq potential depends 
on the heavy quark spin Sg and leads to the spin-spin interaction. The total angular 
momentum is J = j + S^. Thus, for 5-wave mesons there is a doublet of j — 1/2 states 
with y'' = 0^, 1~ (P is the meson parity) which are degenerate in the limit mg — > oo. 
For P-wave mesons there are similarly two doublets of initially degenerate states with 
7 = 1/2 (./^ = 0+, 1+) and j = 3/2 (7^ = 1 + , 2+). The j = 1 /2 levels are expected to be 
broad because they decay in an 5- wave, while the j — 3/2 levels should be narrow since 
they decay in a Z)-wave. 

The heavy-light quark-antiquark potential in configuration space in the limit mg — > oo 
reads as [2] (Vcoui(^) = -(4/3)aj/r, a, = 0.5 for B,D; a, = 0.45 for 8^,0^) 
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The 1 /mq correction to the potential (5) is given by [2] 

= ;^{p[Vcoul(r)+VeonfW]P + ^Coul(^)^ 
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+ 3 [AVcoul(r) + (1 +K)2Ay<f,^(r)] SgS, [. 



Here the prime denotes differentiation with respect to r, L is the orbital momentum, 
and Sg are the spin operators of the light and heavy quarks, S = Sg + Sg is the total spin. 

First Eq. (1) is solved numerically with the complete potential (5) and then the 
corrections (6) is treated perturbatively. The results of the calculation of heavy-light 
meson mass spectra are presented in Tables 3-6 and in Figs. 1-4 [2]. 

From Figs. 1-4 it follows that in the limit mq oo the P-wave doublet with j = 1/2 
lies higher than the one with j — 3/2 (abnormal level ordering), i.e. these doublets are 
inverted. For finite mg the hyperfine splitting of the doublet states makes this picture 



TABLE 3. Mass Spectrum of D Mesons (Ge V). 



State 


Particle Theory [2] 


PDGi2000^ 


CLEO 


DELPHI 






1.875 


1.8693(5) 






15i 


D* 


2.009 


2.0100(5) 






IP2 


05 


2.459 


2.459(4) 






IPi 


Dl 


2.414 


2.4222(18) 


2.425±0.002± 0.002 




IPi 


Dl 


2.501 




2.46llSStl±0.01 ±0.032 


IPo 


D* 


2.438 








2So 


D' 


2.579 












2.629 






2.637(9) ? 


TABLE 4. 


Mass Spectrum of Ds Mesons (GeV). 




Stale 




Panicle 


Theory |2| 


PDG(2000) 


FOCUS 


15o 




Ds 


1.981 


1.9686(6) 








D*s 


2.111 


2.1124(7) 








D*s2 


2.560 


2.5735(17) 


2.5673(13) 


1/^1 




Dsi 


2.515 


2.53535(60) 


2.5351(6) 






Dsi 


2.569 










D*sO 


2.508 










D's 


2.670 






25i 




D*' 


2.716 







much more complicated (some of the levels from different doublets overlap). Only for 
B mesons the pure inversion is preserved. 

In the limit rriQ — > 0° heavy quark symmetry predicts simple relations between the 
spin-averaged masses of B and D states 

- Mdi = Mb,i - Md,i = Mb, - Md, =MB-MD = mb-mc = 3.33 GeV, (7) 

where Mbi = (3Mfij +5Mfi2)/8, Mb = {Mb + 2>Mb*)/A, etc. From Tables 3-6 the fol- 
lowing values (in GeV) for the calculated mass differences can be obtained 

Mbi-Mdi Mb^,-Md,i Mb,-Mds Mb-Md 
3.29 3.30 3.34 3.33 

Thus relations (7) are satisfied with good accuracy. 

The hyperfine mass splittings of the initially degenerate P states 

AMb = Mb2 - Mb^ , Mbi - Mbo ; AMd = Md2 - Md^ , Md^ - Mdq 

should scale with heavy quark masses: AMb = {rnc/mh)AMo and the same for Bs and 
Dg mesons. Our model predictions for these splittings are displayed in Table 7 [2]. 



TABLE 5. Mass Spectrum of B Mesons (Ge V). 



State Part. Theor.[2] PDG(2000) OPAL L3 DELPHI CDF ALEPH 



loo 


R 
1 J 




5 97Qrif'5~> 








B* 


5.324 


5.3250(6) 






IPt 

z 


B% 


5.733 






5.768(8)? 5.732(21) 5.739(13) 




Bi 


5.719 




5.738(9) 


5.71(2) 




Bi 


5.757 






5.670(16)? 


IPo 


B*o 


5.738 




5.839(14) 




2So 


B' 


5.883 








2Si 


1 

B* 


5.898 






5.90(2) ? 






TART V a 


Mass Spectrum of Bg Mesons (GeV). 






State 


Particle 


Theory [2] 


PDG(2000) OPAL 






J- J(j 


Bs 


5.375 


5.3696(24) 








B*s 


5.412 


5.4166(35) 






IPI 


B*s2 


5.844 


5.853(15) 






m 


Bsi 


5.831 








iPi 


Bsi 


5.859 








IPo 


B*sO 


5.841 








2So 


K 


5.971 








2Si 


Bf 


5.984 





CONCLUSIONS 

So one may conclude that for heavy quarkonium mass spectra most of theoretical 
calculations are generally in good agreement with high-precision experimental data. 
This makes possible an accurate determination of fundamental parameters governing the 
dynamics of the heavy QQ interaction. On the other hand, for heavy-light meson mass 
spectra the situation is much more indefinite and unclear in both theory and experiment. 
Not all models predict the P level inversion. The broad j = 1/2 levels are very poorly 
studied due to difficulties in their observation. Much further work is required to remove 
discrepancies between existing experimental data. 
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00-02-17768 and Russian Ministry of Education under Grant No. EOO-3.3-45. 



TABLE 7. Hyperfine Splittings of P Levels (MeV). 



States 


AMd 


mi, ^ 


AMb 


AMd^ 


i^^AMz) 


AMs^ 


IP2-IP1 


45 


14 


14 


45 


14 


13 


IPi - l/'o 


63 


20 


19 


61 


19 


18 
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IS 



mQ— >infimty 



2S1 
2S0 

IPl 
1P2 
IPO 
IPl 



ISl 



ISO 



with 1/mQ corrections 



FIGURE 1. The ordering pattern of D meson states. The mass scale is in GeV. 
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with 1/mQ corrections 



FIGURE 2. The ordering pattern of Dj meson states. The mass scale is in GeV. 
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FIGURE 4. The ordering pattern of Bj meson states. The mass scale is in GeV. 
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